Non-Fermi Liquid behavior in Neutral Bilayer Graphene 
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We calculate the density-density response function and electron self-energy for undoped bilayer 
graphene, within the Random Phase Approximation (RPA). We show that the quasiparticle decay 
rate scales linearly with the quasiparticle energy, and quasiparticle weight vanishes logarithmically 
in the low-energy limit, indicating non-Fermi liquid behavior. This is a consequence of the absence 
of a Fermi surface for neutral bilayer graphene and corresponding larger phase space available for 
scattering processes. Experimental consequences of our results as well as their differences from those 
of single-layer graphene are discussed. 



Introduction — The isolation and subsequent identifi- 
cation of graphene, an atomically thin electron system, 
has led to intense experimental and theoretical intereslii. 
Recent experimental progress^ has also led to techniques 
that enable isolation and study of systems with a small 
number of graphene layers, of particular importance is 
AB-Bernal stacked bilayer graphene^, a system which 
shares some features both with graphene and two 
dimensional electron gas (2DEGs)^, however at the same 
time different from both. Neglecting trigonal warping 
balanced bilayer graphene can be identified as a zero-gap 
semiconductor with quadratic dispersion; for undoped 
bilayer graphene the Fermi energy lies at the neutral 
Fermi point, described as the point where the degenerate 
particle- hole bands meet. Collectively these systems 
can be classified as chiral 2DEGs^. Electron-electron 
interactions in chiral 2DEGs can lead to interesting 
quasiparticle properties, for example quasiparticle ve- 
locity enhancement^ in graphene due to the presence of 
unscreened Coulomb interactions. Most of the physics 
in this paper focuses on the difference in graphene and 
bilayer graphene's chiral 2DEG. 
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In this paper we investigate the quasiparticle proper- 
ties of undoped bilayer graphene due to electron-electron 
interactions. Short-ranged interactions for quadratic 
dispersion in two dimensions are marginal^i^i^ at the 
tree level while Coulomb interactions are relevant; this 
already points to the possibility of non-Fermi liquid be- 
havior. Unlike neutral graphene where electron-electron 
interactions are unscreened. Coulomb interactions in 
bilayer graphene are screened due to the presence of 
a finite density of states, and dynamically generate a 
momentum scale A = me^/e proportional to the inverse 
Thomas- Fermi screening length qtf — 4Alog[4]. Below 
this scale this system resembles one with effective short- 
ranged interactions. Based on the scaling form of the 
density-density response function we demonstrate within 
the Random Phase Approximation (RPA) that there 
is no renormalization to the electron effective mass, 
while the imaginary part of the electron self-energy 
XmT, ~ uj below the screening scale. As a result the 
quasiparticle has a logarithmically vanishing spectral 
weight at low-energies, indicative of non-Fermi liquid 
behavior. This is very different from a 2DEG where 
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FIG. 1: (Color online) Intensity plot of electron spectral func- 
tion at long- wavelengths and low-energies in units of m/X^. 
The solid line corresponds to the noninteracting dispersion, 
m is the electron effective mass and and A — me^/e is pro- 
portional to qtf the Thomas-Fermi screening wave-vector. 
The quasiparticle width (or scattering rate) is proportional 
to its energy and the quasiparticle spectral weight vanishes 
logarithmically at low-energies. See text for details. 



the phase space available for scattering is limited by 
the energy thereby giving energy squared dependence 
for the quasiparticle decay rate, and the quasiparticle 
weight remains finite at the Fermi surface. We argue this 
non-Fermi liquid behavior is a consequence of the ab- 
sence of a Fermi surface for neutral bilayer graphene and 
corresponding larger phase space available for scattering 
processes, which should be a robust result beyond RPA. 
The long-wavelength behavior of the electron spectral 
function A{k,uj) — A{k,—uj) is plotted in Fig. [1] which 
can be compared with ARPES measurements currently 
being employed to study the effects of interactions in 
graphene systems^. 

Bilayer Graphene Effective Model — The low energy 
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properties of Bernal stacked bilayer graphene can be ad- 
equately described by quasiparticles with parabolic dis- 
persion^ exhibiting a Berry phase of 2-k. Neglecting trig- 
onal warping the band hamiltonian for balanced bilayer 
graphene is: 
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where the Pauli matrix acts on the two-degenerate 
(K and K') valleys, k is two-dimensional envelope func- 
tion momentum measured from the two nodal points K 
and K', and are Pauli matrices that act on bi- 
layer graphene's pseudospin (layer) degrees of freedom, 
and a =1, {[) accounts for the spin degrees of free- 
dom. The chirality of bilayer graphene's chiral 2DEG 
is captured by the unit vector h-^ — [cos2Lp^^sai2ipj:) 



where Lpj^ = tan ^(ky/kx). In Eq. ([T]) the field oper- 
ator ^ih^^"^ — ^^'^''f* /A*^'"^ 
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- ,(t>"''' - ,1^"'' - ,(?!>" - ) is a four- 

K+k,t' ^K+k,b' ^K' + k,b' ^K'+k,t' 



component spinor where the low energy sites^ are the top 
(t) and bottom (b) layer sites without a near-neighbor in 
the opposite layer. The effective mass is determined by 
TO = 7i/2u^ ~ 0.054me where v is the single-layer Dirac 
velocity and 71 ~ 0.4eV'^'^*' is the inter-layer hopping am- 
plitude. 

The interaction contribution to the bilayer graphene's 
hamiltonian is layer-dependent: 
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where 



^^■01" ^-02 is the total density per spin. 



^q,a = l/2(/5j.Q — Pq,a) ^hc z-compoucut of the cor- 
responding pseudospin density operator in the K valley, 
with 5q — > —S^ in the K' valley, and v± are the symmet- 
ric and antisymmetric combinations of the interaction po- 
tentials in the same (different) layers Vs — 2Tre'^ /eq{vd — 
Vse~'^'^) with the layer separation d = 0.334nm. 

Response Functions — The finite temperature non- 
interacting response functions can be written as 



in Physically these form factors determine the rel- 
ative weight of the interband (ss' = —1) and intraband 
{ss' = +1) excitation contribution to the density-density 
and pseudospin density response functions. 

For undoped bilayer graphene where the Fermi en- 
ergy lies at the neutrality point (i.e ep — 0), the 
zero-temperature density-density response is completely 
determined by interband excitations where the prod- 
uct ss' = —1. Our low-energy theory has a natural 
high-energy momentum cutoff given by the bandwidth 
kc = v^2to7i, from dimensional analysis of ([3]) it is clear 
that the zero-temperature density-density response func- 
tion Xpp = gm/{2TT)<^pp{q/kc,miVLn/q^), where gm/ [2tt) 
is the constant density of states for parabolic disper- 
sion, and $pp is a dimensionless scaling function. The 
density-density response function is free of any diver- 
gences, and sending the bandwidth kc ^ 00 just gives 
Xpp = gfTi/ {2Tr)^pp{minn/q'^)- Analytically continuing 
iCln —fQ + irj the respective real and imaginary parts of 
scaling function <i>pp areii (y = mfl/q^): 
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where (/ = 4 accounts for the spin and valley degener- 
acy, a = p,S^ denotes the density and pseudospin den- 
sity response functions, s, s' = ± is the chiral band in- 
dex, es{k) = sP /{2m), np{es{k)) is the Fermi-Dirac dis- 
tribution and Oj: - is the angle between the wavevec- 

tors k and k + q. The angular dependent matrix ele- 
ment for the density-density response J-pp {k,k + q) = 
1/2(1-1- ss' cos 20g g_|_-) is different from the pseudospin 

density response J^^^' {k,k + q) = 1/2(1 — ss' cos 20^ ^_|_ -) 
due to the presence of the pseudospin operator 5| oc 



valid for 17 > 0. It is important to note that the 
above expression is a function of a single scaling variable 
y = mil I q^ , contrary to the case of any two-dimensional 
system with a Fermi wavevector kp where it is a function 
of two variables namely mfl/kp and q/kp. This scaling 
behavior is a consequence of the absence of a Fermi sur- 
face for neutral bilayer graphene. 

The imaginary part of the density-density response 
function {TmXpp{q,^) = for > q'^/{4:m)), defines 
the edge of the particle-hole continuum. For interband 
excitations this is just given by the minimum and max- 
imum values of f2 = ^k+q ^~ ^k' minimum en- 
ergy for a particle-hole pair is attained for Oj: ^ = tt, 
this can be seen by completing the square and writing 
min(ri) = l/m[{k — q/2Y + 9^/4]. The excitation spec- 
trum here SI > q'^/{Am) is similar to that of neutral 
graphene Q > v\q\ with the difference coming from switch 
to a parabolic dispersion. 

Within the Random Phase Approximation (RPA)^^: 
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we recover the static screening (fi 
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limit of 



„RPA 



{q, n)) form of Ref. [IJ, with the Thomas-Fermi 
wavevector qrp — g{me'^)/e\og[4:]. Due to the posi- 
tive definite value of 7?,eXpp((7, Jl), bilayer graphene at 
the neutrality point excludes any plasmon excitations. 
The absence of a plasmon mode is not unique for neu- 
tral bilayer graphene and is phenomenologically similar 
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quasiparticle lifetime can be expressed as: 
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FIG. 2: (Color online) Scaling function h+(y) for the imag- 
inary part of the electron self-energy with Coulomb interac- 
tions in the long-wavelength limit plotted as a function of 
y = mu/k^. The inset show the scaling function h-{y) also 
plotted as a function of y = mu/k^ . (see text for details). 
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The efe-dependence can be readily obtained by dimen- 
sional analysis which gives 1 /t^ ^ oc ek', an exact calcula- 
tion yields I/t,,. ^ = 0.1076gTOUgy^efe. The linear energy 
behavior only depends on the fact that $ is a function of a 
single scaling variable y — mQ./q^ , which is a consequence 
of the scale invariance of neutral bilayer graphene, and is 
independent of the detailed behavior of <i>. Bare Coulomb 
interactions would give a different result, however screen- 
ing dynamically generates a new scale qxF, below which 
interactions effectively behave as short ranged. Based on 
the effects of screening and independence of the inverse 
lifetime on the function $ we anticipate non- Fermi liquid 
behavior in the long- wavelength limit. 

The retarded quasiparticle self-energy within RPA can 
be written as: 



to the case of neutral single-layer graphene. In neutral 
single-layer graphene the density-density response func- 
tion vanishes inside the particle-hole continuum thereby 
excluding plasmon excitations, technically different from 
the case of bilayer graphene. Using the continuity equa- 
tion we can relate the optical conductivity to the density- 
density response function. The RPA optical conductivity 
(accounting for spin and valley degeneracy): 



a{VL) = 4 lim ■ 
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(restoring K) is purely real and frequency independent 
similar to the case of single layer graphene. The optical 
conductivity is not influenced by interactions within RPA 
and retains its universal value of e^/2?i. 

Quasiparticle properties at T — — The essential fea- 
ture of a Fermi liquid is encoded in the imaginary part 
of the quasiparticle self-energy (XmS) or inverse quasi- 
particle lifetime, which for a homogeneous 2DEG gives 
ImE ~ {e — eF}'^ log |e— epl- In the case of neutral bilayer 
graphene the absence of a Fermi surface coupled with 
the effect of the interband excitations invalidates such a 
description. In this section we calculate the imaginary 
part of the quasiparticle self-energy for neutral bilayer 
graphene and show that XraE oc for — )■ indicating 
non- Fermi liquid behavior. 

First let us examine Fermi's golden rule expression for 
short-ranged (g-independent) interaction Ueff. This is a 
rather crude approximation but as we argue below still 
gives the essential features of the predicted non-Fermi 
liquid behavior. Based on Fermi's golden rule the inverse 
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following the the line and residue decomposition of Quinn 
and Ferrell^^. The quasiparticle self-energy within RPA 
remains diagonal in the particle-hole basis. It can be 
shown that the line contribution is purely real and does 
not contribute to the imaginary part of ImEIjf* which 
for {uj>0): 
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where £(g, fi) = 1 + w+((7)xpp(g, fi). In the above ex- 
pression we have neglected the contribution of V- as it is 
logarithmically suppressed once screening effects are ac- 
counted for within RPA. Dimensional analysis of ([9]) im- 
plies that ImY,+ {k,u) = \k\^f{mu)/\k\^,k/qTF)^ where 
/ is a two variable function. Interactions introduce an 
inverse length scale A = me^ je which turns out to be the 
same order of magnitude as the bandwidth cutoff fee- In 
the long- wavelength limit we find that: 



ImTrf{k ^ 0,w ^ 0) 
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(10) 

The scaling function h-f-(muj/k'^) was numerically at- 
tained and is plotted in Fig.[2l The theta function in (fTO]) 
comes from the particle-hole continuum and is indepen- 
dent of the nature of interactions. Using the symmetry 
relations ImE^'^*(fc, —uj) = ImY7'/*-{k^ lu) where s = — s it 
is clear to see that ImE!j5*(fc, —lu) gives a similar expres- 
sions as PU)) with a different scaling function h_{muj /k"^) 
plotted in the inset of Fig. [21 The residue contribution to 
the real part of the retarded electron self-energy TZeT:^* 



[6(0; -e(g))], (9) 
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yields a similar expression, however the contribution due 
to E^™*^ is more singular thereby dominating in the long- 
wavelength limit. We find that (w > 0): 

neJ:r{k,u;) = .^(log [^])2-ii^(log + 

(11) 

where represent the subleading terms. The expres- 
sion for Lu < can be attained by exploiting the sym- 
metry relations 7^eS7*(fc, -tj) -7^eI;^'^*(fc, w). The 
quasiparticle spectral weight: 
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vanishes logarithmically, whereas the effective mass m* 



mi 



1 -a^7^eS!;'=*(fc,w) 



1, 



(13) 



remains finite and is not renormalized by interactions. 
The long-wavelength behavior of the spectral function 
A{k,u!) = A+{k,uj) + A^{k,u!) plotted in Fig. [1] with: 



\lo - Ek - nej:i''*{k,Lj)\^ + 



|XmI]^'=*(/c,cj)|2 
(14) 

was calculated from the leading order behavior of the 
electron self-energy and neglecting regular contributions. 
Symmetry relations for the electron self-energy dictate 
that A{k,uj) = A{k,-uj). 

The linear dependence of the JmS on the quasipar- 
ticle energy predicted above is different form the case 
of neutral single layer graphene. In neutral single layer 
graphene due to the lack of screening associated with the 
Dirac point the Fermi velocity develops a logarithmic 
enhancement^. For neutral graphene this logarithimic 
velocity enhancement implies that ImY^iuj) ^ oj / {loguj)'^ 
which is smaller than Im'E{uj) lo. In contrast to 
single layer graphene interactions in bilayer graphene are 
screened, as we have shown and has been pointed out 
in the literaturei^^ with Thomas-Fermi screening, that 
the quasiparticle dispersion in bilayer remains quadratic. 
Most of our analysis of quasiparticle properties in 
neutral bilayer graphene has relied on the fact that 
Xpp(g, 5^) oc Ld^ (i.e. it has a scaling behavior after 



sending the bandwidth fcc oo). In the next section 
we identify this regime of non-Fermi liquid behavior at 
finite temperature. 

Regime of Non-Fermi liquid behavior — Weak inter- 
layer hopping leads to trigonal warping of the band struc- 
ture in bilayer graphene-. The temperature associated 
with this effect can be estimated by calculating the en- 
ergy scale at which trigonal warping effects compete with 
the quadratic dispersion kept within our model. Using 
the bare parameters of graphite 73 ~ 0.l7cr^ we can 
estimate that the temperature below which the trigo- 
nal warping effect becomes relevant is Ti ~ 40 K. Be- 
low this scale the electron dispersion crosses over from 
quadratic to linear, and the system behaves like single- 
layer graphene, and our results no longer apply. Recently 
there has been discussions of interaction-driven sponta- 
neous symmetry breaking in neutral bilayer graphene in 
the absence of trigonal warping, due to the marginal rele- 
vance of weak short-range repulsive interactions^!^. Since 
the interaction is only marginally relevant, the transi- 
tion temperature Tc into the possible broken symmetry 
phases are exponentially small. The non- Fermi liquid be- 
havior discussed here thus applies to temperatures above 
the higher of Tc and Ti. Our analysis can be extended 
to finite temperatures, again dimensional analysis dic- 
tates that the polarization function have the scaling form 
^{mVl/kpT^mQ./q^). From a simple scaling analysis of 
the electron self-energy we conjecture that for tempera- 
tures T > max{Ti,Tc): 
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UJ > ksT 
ksT > LU 



(15) 



indicative of non-Fermi liquid behavior. In this paper 
we have analyzed the regime k < qxp, in the opposite 
regime k > qxp one can show that the non- Fermi behav- 
ior becomes even more pronounced as screening effect is 
less significant. 
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